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Abstract -We consider the rapid quench of a one-dimensional strongly correlated supersolid 
to a localized density wave (checkerboard) phase, and calculate the first-order coherence signal 
following the quench. It is shown that unique coherence oscillations between the even and odd 
sublattice sites of the checkerboard are created by the quench, which are absent when the initial 
state is described by a Gutzwiller product state. This is a striking manifestation of the versatility 
of the far-from-equilbrium and nonperturbative collapse and revival phenomenon as a microscope 
for quantum correlations in complex many-body states. For the present example, this opens 
up the possibility to discriminate experimentally between mean-field and many-body origins of 
supersolidity. 



The single-particle density matrix, from which the first- 
order coherence of the system is defined, represents a sensi- 
tive probe of the correlation properties of interacting quan- 
tum many-body systems. This becomes especially appar- 
ent in far-out-of-equilibrium dynamics when one rapidly 
(nonadiabatically) quenches from one quantum phase to 
another, which has recently led to intensive research into 
the properties of the resulting dynamical quantum phase 
transitions PQ. By now the most familiar nonequilibrium 
phenomenon occurring in such transitions are the collapse 
and revival oscillations of first-order coherence, created 
after sweeping from the supcrnuid to the Mott phase in 
optical lattices. When only conventional contact interac- 
tions are present, this has been first observed in the classic 
paper [2j . Experiments on the superfluid-Mott transition 
are now entering a high-precision age [3] , enabling the de- 
tection of many-body correlation functions and their tem- 
poral behavior on the single-particle and single-site accu- 
racy level. Stimulated inter alia by this fact, various facets 
of the collapse and revival phenomenon were theoretically 
investigated recently in a number of papers, e.g. in [4]. 

The extended Bose-Hubbard model, realizable experi- 
mentally via cold atoms or molecules with large dipolc mo- 
ments stored in optical lattices [S], leads to a significant 
extension of the phase diagram of bosons on an optical 
lattice, including supersolid phases. By definition, for su- 



pcrsolids, the order parameter amplitude, which is homo- 
geneous for a supcrfluid ground state, becomes spatially 
modulated, caused in the Bose-Hubbard model by two- 
body interactions coupling neighboring lattice sites with 
each other. While the implementation of supersolids in 
optical lattices cannot answer the intensely debated ques- 
tion whether they occur in real solid-state systems like 
4 Hc, where the existence of (growth-induced) defects in 
the crystal complicates the interpretation of experiments 
6 10], the microscopic and thus many-body origin of the 
supersolid phase can be understood better in such an ar- 
tifically engineered system. 

In the following, we derive a unique signature for the im- 
print of the initial many-body correlations contained in the 
supersolid state from the far-out-of-equilibrium collapse 
and revival oscillations of the extended Bose-Hubbard 
model after a quench. The quench envisaged is the ana- 
log of the superfluid-Mott quench in that a delocalizcd 
phase-coherent (here supersolid) phase is quenched to an 
incoherent and localized, here density wave (DW) phase. 
We shall demonstrate that the final time-dependent co- 
herence signal can be used as a sensitive probe of the ini- 
tial strongly correlated supersolid state and therefore of a 
possible many-body origin of supersolidity. We show that 
even for very small density modulations of the initial su- 
persolid ground state, which might be difficult to detect 



p-1 



directly in experiment |lllll2j , a distinct first-order coher- 
ence signal appears after the quench, providing a unique 
signature of supersolidity caused by both nonlocal inter- 
action couplings and strong correlations. 

We consider here that quantum fluctuations are strong 
in the initial ground state, for a quasi-one-dimensional 
(quasi-lD) Bose gas in a ID optical lattice directed along 
the axis [T3JQ3]. We determine the ground state by ex- 
act diagonalization (ED), and contrast the correlation re- 
sponse after the quench thus obtained with the temporal 
correlations when the initial ground state is assumed to be 
describable by a Gutzwiller product state (GW) [5]. It is 
stressed that we propagate both ED and GW initial wave- 
functions exactly. We note in this regard that applying the 
time- dependent GW mean-field approach, used to describe 
small oscillations on top of the ground state, cf. [TS], to 
situations far from the (instantaneous) ground state (e.g. 
for collapse and revival oscillations), leads to both quanti- 
tatively and qualitatively incorrect predictions in the ex- 
tended Bose-Hubbard model [IB]. By exact propagation, 
we conclusively demonstrate the strong dependence of the 
dynamical correlation response after the quench on the 
presence of many-body correlations in the initial state. 
The final time-dependent correlation functions can there- 
fore serve as a sensitive probe of whether the initial state 
correlations are correctly captured by a specific ground- 
state wavefunction ansatz for the supersolid state. 

The extended Bose-Hubbard Hamiltonian, describing 
bosons with dipolar interactions in the lowest Bloch band 
of an optical lattice, contains the nearest neighbor tunnel- 
ing rate J, onsite interaction coupling U, and finally the 
nearest neighbor (NN) coupling V. Assuming a homoge- 
neous system (no external trapping potential) , it reads 

M M y M 

H = -J 4 & j + ~2 ~ + ~2 Aifl i' 

<ij> i <ij> 

(1) 

where M denotes the number of sites and < ij > indicates 
a summation over NN (we will provide below a discussion 
of the influence of including off-site interaction couplings 
beyond NN) . The onsite and offsite couplings U and V are 
proportional to the dipole moments squared and depend 
on integrals over products of the Wannier orbitals for the 
lowest Bloch band. A detailed derivation of the extended 
Bose-Hubbard Hamiltonian for dipolcs in an optical lattice 
can be found, for example, in [5]. We assume the dipolcs 
to be oriented perpendicular to the axial direction such 
that the NN coupling is repulsive, V > 0. We use U = 1 
as unit of the energy and 2-k /U as the unit of time in what 
follows (H = 1); the off-site coupling is fixed at V = 0.9; 
this choice of V < U is motivated by the fact that then 
small population imbalance supersolids can be obtained at 
relative tunneling rates J /(nil) <C 1, also cf. (T9"] . 

Previous studies into the extended Bose-Hubbard 
model, covering the whole range of fillings and couplings, 
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Fig. 1: Top: Artist's view of the supersolid ground state in 
a laterally tightly confined quasi- ID gas with axially superim- 
posed optical lattice (solid line). For representation purposes, 
the density imbalance Iss is chosen large. Bottom: Correla- 
tion function g\i/n of initial state for M — 12 and N = 8 from 
ED with n = 2/3 and J = 0.05 (open circles); solid lines are a 
guide to the eye. 



have found various other phases like superfluids and lo- 
calized Mott and DW phases besides supersolids [20l[2T] . 
We concentrate in what follows on supersolid initial states 
and final (after the quench) DW states. 

The ground-state density modulation in the supersolid 
state for the lattice model ((TJ) is related to a structural in- 
stability, developing when the excitation spectrum above 
the superfluid ground state, at the transition to the super- 
solid state, touches the quasi- momentum axis at q — ir/a 

P21CE71CEB]. This leads t0 " A " ( with indices i = A = 
1, 3, 5, . . .) and "B" (with indices i = B = 2, 4, 6, . . .) sub- 
lattices following a checkerboard pattern associated with 
the density modulation of the supersolid ground state at 
wavevector q = ir/a. A single scalar measure distinguish- 
ing the supersolid state from the superfluid state is there- 
fore the population imbalance 



\nA - n B \ 



(2) 



where Tia and tib are average fillings of sublattice sites; 
hence Iss = for a superfluid and Igs ^ for a supersolid 
state. 



The phase and modulus of the single-particle density 
matrix 



(7^; = {a k ai) = cxp[i6> 



(3) 



embody density and phase correlations between sites, 



and lead to first-order correlations gij 
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|(ajaj}| cos[0jj]. We calculated the ground-state correla- 
tion function corresponding to the Hamiltonian ((T|) for a 
supersolid state in a lattice with M = 12 sites and N = 8 
respectively N = 16 particles, i.e., at fillings n = 2/3 
and h = 4/3. For this system size, the determination of 
both initial ground state and temporal evolution can be 
performed by ED routines, employing periodic boundary 
conditions. In addition, we have verified the qualitative 
robustness of our results by comparing with sample cal- 
culations for M = 10 sites at fillings n — 4/5 as well as 
M = 14 sites at filling n = 9/14 for ground state and 
quench dynamics, and M = 16 sites at n = 5/8 for the 
supersolid ground state. 

There is the inherent lattice symmetry leading to degen- 
erate SS states which have either ua > ng or tla < ng. 
We choose one of these states, say with ua > ng, by the 
following procedure. We perform a large sequence of nu- 
merical calculations with different initial seeds for many- 
body states in the ED calculation, to obtain a statistical 
distribution of various values of Iss with different degen- 
erate states. From this distribution, we choose the initial 
state as the one with the most probable value of Iss which 
fulfills tia > ng. 

As a result, we obtain for the supersolid ground state 
that, at constant phase 8jk, first-order correlations are 
strictly positive, decaying with distance in the power-law 
fashion characteristic of low-density quasi- ID Bose gases 

na, cf.Fig.ni 

Next we let the tunneling abruptly go to zero in a quasi- 
instantaneous quench, J — » 0, which would, in an adia- 
batic transition, transfer the supersolid to the localized 
and insulating DW phase. The main feature we observe 
is that there are two different coherence signals obtained 
after the quench, depending on the distance between sites 
considered being \i—j\ = even or odd. This is a first crucial 
difference from the uniform (long-range) collapse and re- 
vival in the superfluid-Mott quench OS]. The "AA" (and 
equivalently "BB" ) correlations between the sites belong- 
ing to the same sublattice (\i — j\ =cven) show (partial) 
collapse and revival oscillations, i.e. with strongly reduced 
amplitude, at constant phase 9aa (Obb), cf. Fig.[2][b). On 
the other hand, the "AB" correlations belonging to sites of 
different sublattices — j\ =odd) show periodic changes 
from positive to negative values of correlations, see Fig. [21 
while the average populations Tla are n b are both constant 
in time (implying that the phase difference switching does 
not lead to a redistribution of occupation numbers) . The 
phase difference Bab in the latter case thus switches from 
values close to zero to close to n, for prolonged time peri- 
ods of order the ordinary revival time scale 2n/U [we are 
using a representation of the phase difference which takes 
into account that it is defined modulo 2ir, i.e. in particu- 
lar identify 2ir and 0] . The first order correlation function 
9AB{t) goes through zero when #ab switches from zero to 
7r and vice versa. For these out-of-cquilibrium states far 
from the DW ground state, we coin the term ir-states [23] . 

We stress that for an ordinary superfluid to Mott quench 
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Fig. 2: Coherence and phase-difference signals obtained in ED 
(we fix i = l,j = 4 in <?ab and i = l,j = 5 in Qaa)- Left col- 
umn: gAs/n (a) and gAA/n (b) for filling n — 2/3 at J = 0.05 
(i.e. with initial correlations shown in Fig.QJ. Right column: 
Corresponding phase difference signals, (c) shows the gAB and 
Oab results for higher filling n = 4/3 and at J = 0.08 (note 
the different time axis). Imbalance Iss = 0.0136 in (a),(b) and 
Iss = 0.009 in (c). 

at V — 0, no 7r-state is obtained, even in the present 
strongly correlated low-filling ID lattice case. Hence a 
nonvanishing off-site coupling V in the Bose-Hubbard 
Hamiltonian ([1]) is crucial for obtaining the large nega- 
tive AB correlations characteristic of the 7r-state. Further- 
more, we also observed that when the initial imbalance Iss 
becomes large (achieved when we take V > U at small J), 
and therefore the initial supersolid state is less strongly 
correlated across different sites, that the tendency to form 
a 7r-state disappears. In the extreme case of Iss — * 1 
(one sublattice empty), the phase simply changes linearly, 
9 A B °c t. 

We now contrast the occurrence of 7r-states, obtained 
starting from a strongly correlated initial state, with the 
first-order correlation signal from an initial supersolid 
state assumed to be given in GW mean-field theory. The 
many-body state can then be written as the product state 

m 

I GW > = I1 E fA{nA)fB(n B )\n A ,n B ), (4) 

A,B UA,n B 

where J\ A B is a shorthand for taking the product over 
all sites, and the indices A and B run over all M/2 cor- 
responding sites of either sublattice. The sublattices are 
separated here explicitly, with j a and /b in general being 
different distribution functions for cither sublattice. 

When the initial state is given in the product form (j4j, 
the exact time evolution of AB correlations in the limit 
of zero tunneling, J = 0, after the quench can be written 
in the form of Eq. ([5]) (with a similar expression for the 
AA and BB correlations) , where NN I designates all sites / 
which are nearest neighbors to the sites with given indices 
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{a\a B )= fA(n A )fB(n B )fX(n A + l)f* B {n B - l)y/n B {n A + iy tu ( n *- n °) 

n A ,n B 

x nEi^)i 2 i/w)i 2e ' w(E[ " ]n '" E[B,ini) - ( 5 ) 

NN! ni 



A and B, and we assume for simplicity that the distance 
\A - B\ > 2 US]. The notation [Al] and [Bl] indicates 
sums over all NNs I at the sites A and B, respectively, 
and f(A,l) means f A (ni) where I is a NN to A and for 
f(B,i) correspondingly. Wc therefore see that in case the 
initial state can be written in Gutzwiller product form (QJ , 
the final-state time-dependent correlations solely depend 
on the on-site initial number distribution functions in the 
A,B sites, and in particular on their widths and centers. 
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Fig. 3: Coherence signal gAB(t) between sublattices obtained 
when propagating exactly a GW supersolid initial state after 
a quench to J = 0, for the indicated values of filling and tun- 
neling [cf. for the ED results in Fig. [2] (a) and (c), left col- 
umn]. Top left has the same Iss = 0.0136 like the ED result of 
Fig-UJa), while on the right Iss = 0.5339. Bottom for larger 
filling, left again Iss = 0.0136 and right Iss = 0.25. 



We time evolve the initial GW state ((3]), found upon 
minimizing (GW| H — fxN |GW), where \i is the chemical 
potential, using the exact formula (O , and show the results 
in Fig. [3] While there is some negative amplitude of g AB 
in the large-imbalance low-filling case (cf. the right upper 
plot in Fig. [3]), the highly symmetrical 7r-state of Fig.[2](a) 
and (c) is not obtained. For illustration purposes, we dis- 
play in Fig.[3jthe evolution results for the same low fillings 
like in Fig. [2] for which the GW approach naturally cannot 
be expected to give the correct ground state, in particu- 
lar for small population imbalance. Generally speaking, 
a much larger J is necessary to obtain in GW the same 
Iss as compared to ED at the same filling, which empha- 
sizes the strong correlation of a low-filling supersolid in 
ID. We have furthermore confirmed that the GW correla- 
tion response is predominantly determined by the imbal- 
ance Iss, and at the same values of the latter quantity 
remains qualitatively similar at larger filling n. We also 
note that the maximally possible sublattice imbalance for 
the supersolid ground state decreases with increasing n, 



so that at larger filling, where the non-number-conserving 
GW theory becomes asymptotically exact, the tendency 
of developing negative <?ab as shown in the right row of 
Fig.[3j vanishes completely. The fact that with an initial 
GW product state large negative AB correlations are not 
obtained serves to illustrate that both large nonlocal in- 
teraction coupling and strong correlations are necessary 
to obtain the large negative values of g AB shown in Fig. [5] 
and thus an (almost) perfect 7r-state. 

We also investigated in addition whether the 7r-state is 
robust under increasing the number of off-site interaction 
couplings taken in the extended Bose-Hubbard Hamilto- 
nian ([l]). For a dipolar interaction, the off-site interaction 
Hamiltonian for polarized dipoles reads 



H. 



dd 



M A 
V sr-^ UiTlj 



■3 



(6) 



We are truncating the coupling summation at the distance 
\i — j\ = A, where the NN off-site coupling in (Q]) corre- 
sponds to A = 1. The ED results for the time dependent 
9 AB , for three different cutoff lengths, and the identical 
parameter set of Fig. [5] (a), are shown in Fig.SJ We see 
that the 7r-state, consisting in the switching of the relative 
AB phase between and n, still occurs, although with 
a different period. The period first increases when going 
from A = 1 to A = 2, and then decreases again for A = 3. 
This behavior is due to the ABAB sublattice structure 
of the supersolid initial state. The 7r-state period thus 
encodes nontrivial information on the interaction range 
chosen for the extended Bose-Hubbard model in a quench 
with a periodically modulated initial state. 
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Fig. 4: 7r-state switching of the relative sublattice phase 9ab 
for different cutoff lengths A of a dipolar off-site interaction. 
J = 0.05, n = 2/3 (M = 12, N = 8); Solid black line A = 1, 
dashed red line A = 2, dotted blue line A = 3. 
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The 7r-states can be detected experimentally, for ex- 
ample, by the time-dependence of peaks in the quasi- 
momentum distribution function 

^(*)4E(^i)^ (M . (7) 
k,l 

which correspond to the visibility of interference peaks in 
the real-space density detected by time-of-fiight experi- 
ments [TT]. Pronounced peaks are located at q = and 
q = n/a, cf. Fig.O whose time dependence should be eas- 
ily discernible. This is due to the fact that gAA(t) [<?bb(£)] 
and gAB{t) strongly differ on average and have sufficient 
oscillation amplitude even for small population imbalance 
Iss, cf- Fig. [21 which leads to significant temporal vari- 
ation and amplitude of n q (t). As shown in Fig.JSJ the 
appearance of the 7r-state leads to a sizable increase of the 
peak at q = n/a, while it decreases the q = peak. 
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Fig. 5: Quasi-momentum distribution function n q (t) for two 
different times. The time t' is defined by the influence of the 
7r-state on the quasi-momentum distribution being most sig- 
nificant, corresponding to the negative minima of the correla- 
tion function Qab- Parameters are chosen identical to those of 
Fig.[l(a). 



We have presented a study of ID supersolids on opti- 
cal lattices undergoing a quench to a density wave phase, 
and observed that a final nonequilibrium 7r-state signa- 
ture develops, characteristic of strong correlations in the 
initial state. By contrast, the 7r-state does not occur when 
the inital state is described by a Gutzwiller product state 
ansatz. This is a prime example of the high sensitivity of 
far-out-of-cquilibrium many-body dynamics on initial cor- 
relations. As a corollary, even when a Gutzwiller product 
state can reasonably well emulate certain properties of the 
equilibrium ground state, it can potentially fail tremen- 
dously as an initial state for many-body evolution. Cor- 
relation functions after a quench can therefore serve as 
a "microscope" of the initial correlations, in which the 
far-out-of-equilibrium conditions magnify the many-body 
character of a state or the lack thereof. In the present 



case, this offers the possibility to detect a strongly corre- 
lated supersolid state of matter even at very small density 
modulations, and thus to determine accurately the phase 
transition point to the supcrfluid state in which no 7r-state 
occurs. 
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